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We present a detailed study of the Nernst effect N^x in (TMTSF)2PF6 as a function of tempera- 
ture, magnetic field magnitude and direction and pressure. As previously reported there is a large 
resonant-like structure as the magnetic field is rotated through crystallographic directions, the Lebed 
Magic Angles. These Nernst effect resonances strongly suggest that the transport of the system is 
effectively "coherent" only in crystallographic planes along or close to the applied field direction. 
We also present analytical and numerical calculations of the conductivity and thermoelectric tensors 
for (TMTSF)2PF6, based on a Boltzmann transport model within the semi-classical approximation. 
The Boltzmann transport calculation fails to describe the experiment data. We suggest that the 
answer may lie in field induced decoupling of the strongly correlated chains. 



PACS numbers: 74.70.Kn, 72.15Nj, 72.15.Gd 



I. INTRODUCTION 

(TMTSF)2PF6 is a quasi-one-dimensional electronic 
system, which displays various ground states rang- 
ing from triplet superconductoriS*^ to spin-density 
wave(SDW) insulator, depending on pressure, temper- 
ature and magnetic field (for a review, please see Ref0). 
(TMTSF)2PF6 consists of plate-like TMTSF molecules 
which stack with strong wavefunction overlap in chains. 
The intrachain bandwidth is ~ leV while the inter- 
chain couplings give anisotropic bandwidths of O.leV and 
0.003ev in the approximately orthogonal directions. In 
the 'metallic' phase under moderate magnetic field, a 
fascinating phenomenon, the so-called Lebed Magic An- 
gle Effect (MAE) was discovered^iSiZ, after Lebed's initial 
predictior&S. The first manifestations of these MAEs 
were sharp resistance dips when the magnetic field was 
aligned at inter-chain directions in real space (lattice 
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In reciprocal space a field along the magic 
angles induces electron motion along commensurate k 
space orbitsii. Despite many theoretical efforts to de- 
scribe the magic angle efFectsiSiiiiiSiiiiiiiSiiSiil, there is 
as yet no satisfactory explanation. Most of the theories 
focus on the semiclassical motion of electron on the open 
Fermi surface derived from single particle band structure. 

Recently, a giant Nernst effect was discovered in 
(TMTSF)2PF6. As the magnetic field was rotated to- 
ward a magic angle the Nernst signal increased then de- 
creased toward zero, changed sign at the magic angle and 
continued in an inverse manner. The result is a sharp 
resonant-like structure.--^ The magnitude of the Nernst 
signal at IK is at least 3 orders larger than what we ex- 
pected from simple (Drude) estimates. The sign change 
of the Nernst effect at the magic angles strongly sug- 
gests that the transport involved in the Nernst effect is 
effectively 2-dimensional at these commensurate angles. 
Both the sign change at the magic angles and magni- 
tude of the signal are not yet explained, but the effect 
appears generic for these materials. The giant resonant 
Nernst voltage has recently been observed in the sister 
compound (TMTSF)2C104ii^ Present phenomenological 



models for the sign change involve field induced inter- 
plane decouplingiiS Although there is some experimental 
evidence for this decoupling, there are not yet theoretical 
models which rigorously demonstrate this phenomenon. 
Giant Nernst signals have also been seen in high transi- 
tion temperature superconducting (HTc) cuprates where 
a model invoking superconducting vortices and 2D super- 
conducting phase coherence has been successful. 
A similar model has been proposed for (TMTSF)2PF5 — . 
However, to apply this idea in (TMTSF)2PF6 is quite 
controversial. On one hand it naturally explains the 
large Nernst signal with undetectable thermopower sig- 
nal, predicts a particular sign of the Nernst effect con- 
firmed by experiments and qualitatively explains some 
aspects of experiments. On the other hand, it predicts 
a large superconducting fluctuation region in the phase 
diagram, which is absent in other measurements. Most 
of the superconducting properties in (TMTSF)2PF6 have 
been understood within a conventional mean field BCS 
picture. We will explore the possibility of this vortex 
Nernst effect with more experimental detail in a subse- 
quent paper. 

What sorely hampers progress in understanding these 
unusual magic angle phenomena are the lack of measure- 
ments other than charge transport. The Field Induced 
SDW (FISDW) and MAE have been observed primarily 
in the charge channel by transport measurements. The 
FISDW has been more thoroughly explored with magne- 
tization, magnetocaloric effect and spin relaxation stud- 
ies. Magnetic torque measurements on (TMTSF)2C104 
suggest there is a thermodynamic component to the 
MAE^. While a thermodynamic probe is an obvious 
choice for establishing the presence of unknown phases 
or fluctuations in (TMTSF)2PF6, the high pressure envi- 
ronment makes a measurement of specific heat or dc mag- 
netic susceptibility impractical. Recently, the ''''Se NMR 
spin-lattice relaxation rate measurements^^ at different 
magnetic field orientation show no evidence for either a 
spin gap or a single particle gap. Furthermore, there is 
no evidence for an enhancement of the FISDW transition 
temperature. This strongly suggests that neither FISDW 
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ordering nor fluctuations are likely to be responsible for 
the MAE. The dramatic contrast between the charge 
channel and the spin channel at MAs suggests that spin 
and charge degrees of freedom may be decoupledi ^'^i^^ 
The thermodynamic and suggested coherent-incoherent 
transitions would therefore be the result of interaction 
and correlation effects due to subtle changes in the elec- 
tronic wavefunctions and density wave susceptibilities. 

Before speculating further on exotic mechanisms for 
the giant Nernst resonances and other MAEs in transport 
it is necessary to see what conventional transport theory 
will yield. Although Boltzmann transport calculations 
as a function of magnetic field magnitude and direction 
have been performed for resistance^ ^ there has been no 
such study for the thermoelectric transport coefficients. 
Such calculations are one of the main contributions of 
this paper. 

We divide our presentation into two sections. The first 
section focuses on the Nernst experiments. We present a 
detailed study of the Nernst effect N^^in (TMTSF)2PF6 
at various pressures, magnetic fields and temperatures. 
The second part presents both numerical and analytic 
calculations of Boltzmann transport in the relaxation 
time approximation with realistic band parameters . We 
then compare the calculations with our experiment data. 



II. NERNST MEASUREMENT 
A. Method 

Figure n shows the experimental setup for the Nernst 
effect Nza; measurement with the temperature gradient 
along the a-axis and voltage measured along c-axis. 3 
pairs of Au wires were attached to the opposite ab-planes 
of the sample by silver paint for both resistance measure- 
ments(the end pairs) and the Nernst measurements (the 
middle pair). The Au wires were attached to the al- 
loy wires (Phospher Bronze) fed through the pressure 
cell base. Here we used low thermal conductivity al- 
loy wires instead of Cu wires to minimize the possible 
transverse temperature gradients. A miniature heater 
was placed on top of the sample to establish a small tem- 
perature gradient along the a-axis. Two thin film RuO 
thermometers were used to measure the temperature dif- 
ference. The thermoelectric voltage is measured by a 
Keithley 182 Nanovoltmeter. The heater was turned on 
and off for a few cycles for signal averaging, A linear- 
fit-extrapolation method was used to accommodate the 
slow drift of the baseline signal. The resistance was 
measured by a conventional 4-probe low frequency lock- 
in technique. The magneto-resistance R^^and Nernst sig- 
nals N^jjwere measured simultaneously. 

Fig. 121 shows a typical angle dependence of the Nernst 
signal N^^in (TMTSF)2PF6, obtained at 1 K, 6 Tesla and 
13 kbar. The magnetic field was rotated from -40° to 50° 
with respect to the c*-axis.(See Ref 18 for definitions of 
c*, c' etc.) The maximum Nernst signal is about 100 




FIG. 1: The measurement setup of Nernst effect ^zx'- Three 
pairs of Au wires were attached to the sample along the a axis 
on the opposite sides of the ab-planes of the sample. Two 
RuO thin film resistance thermometers were placed next to 
both ends of the sample to measure the temperature gradient 
generated by a miniature heater on the top. The middle pair 
of leads was used for Nernst voltage AVz pickup. The other 
two pairs of leads were used for 4-probe interplane (c-axis) 
resistance measurement. The magneto- resistance, Rzz, was 
measured simultaneously with the Nernst effect. 



HV/K, found at approximately 3° - 4° off c' {0^,, = 7°). 
The Nernst coefficient is of the order of 10 /iV/K-T. As 
far as we know, this value is much larger than the Nernst 
effect observed in any other metal. The angular depen- 
dence of the Nernst signal agrees well with our previous 
thermoelectric measurement in a different geometry. -'^^ To 
study the temperature, field and pressure dependence of 
the giant Nernst effect, we fixed the magnetic field ori- 
entation at 3° off a magic angle(c' or -IL) . First, let's 
discuss the sign of the Nernst effect in (TMTSF)2PF6. 
This is very important for the vortex Nernst model. 
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FIG. 2; Angle dependence of the Nernst signal Nzx and c- 
axis resistance Rzz were measured simultaneous at 1 K, 6 
Tesla and 13 kbar. The thin line data are Rzzifi)- The open 
circles "o" data are the Nernst data Nzxifi)- The solid line 
is a guide to eyes. The Nernst resonances are well aligned 
with the magic angles marked by the resistance dips. Here 
±1 correspond to inter-chain directions c'±b'. 
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B. The sign of the Nernst effect 

In the vortex liquid phase of a type II superconductors 
vortices flow down the temperature gradient,v|| (— VT) 
and generate an electric field E = B x v transverse to 
the temperature gradient — VT according to Josephson 
relationSSiSS. Therefore, the sign of the vortex Nernst 
effect is fixed by VT x B. In general, the Nernst ef- 
fect of an electronic system can have either sign depend- 
ing on details of the band structure. To determine the 
sign of the Nernst signal, we noted the orientation of the 
sample and leads, and placed an alignment mark on the 
base/feedthru of the pressure cell and on the cell body. 
We assume the alignment mark doesn't change much on 
pressurization. To get the Nernst sign correct we need 
to know the orientation to better than 90°. We see the 
magic angles where we expect them to be to ^ 15° We 
observed no orientation variation when the pressure was 
increased in the same pressure cell. Our measurements 
show that the sign of the Nernst effect is consistent with 
the vortex Nernst model (but certainly does not prove 
it). 



C. Temperature dependence 
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FIG. 3: Simultaneous measured temperature dependences of 
the c-axis resistance R^^ (thin lines) and the Nernst signal 
Nza; (filled symbols) at 8 kbar for 3° off c'. The Nernst signals 
rise to a maximum around 1 K, then fall exponentially and are 
unmeasurable below ~200 mK. The decrease of the Nernst 
signal is correlated with the upturn of the resistance as T 
decreases. The upturn of R^z indicates the Metal-FISDW 
phase transition. 

Fig|2|shows the temperature dependence of the Nernst 
effect at c' at different field values: 4, 6 and 7.5 
Tesla for 8 kbar pressure. The Nernst signals rise grad- 
ually to a maximum around 1 K as the temperature de- 
creases, then fall off roughly exponentially at lower tem- 
perature and are unmeasurable below ^ 200 mK. Clearly, 
the Nernst signal is non-linear with magnetic field be- 
low 2 K. The temperature dependence of the Nernst sig- 



nal at 7.5 Tesla agrees with the previous Nernst data 
from the thermopower measurement.^® At high magnetic 
fields, (TMTSF)2PF6 enters the FISDW insuhing phase 
at a critical temperature Tc(II) determined by the sharp 
rise of resistance R^j; measured simultaneously. In the 
"metallic" phase, the Nernst effect at c' and ±1 Lebed 
angles have similar temperature dependence, except the 
magnitude is much larger for c'. However, the presence 
of the FISDW phase seems to suppress the Nernst signal 
at c'. As shown in Fig|31 for 7.5 Tesla data, a sudden 
decrease of the Nernst voltage happens at the onset of 
FISDW transition. 
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FIG. 4: Main panel: Temperature dependence of the Nernst 
signal Nza^measured at 3°away from -1 Lebed angle at 8 kbar. 
The magnetic fields are: 7.5, 7, 6.5 and 6 Tesla respectively. 
The arrows mark the FISDW transition temperature obtained 
from resistance measurements Rzz(T) (not shown for clar- 
ity) at various fields. Insert: Temperature dependence of the 
Nernst signal N^ajand resistance RzzOf -1 Lebed angle at 13 
kbar for 7.5 Tesla magnetic field. The dash line marks the 
FISDW transition (~ 350mK). An enhancement of Nernst 
voltage was found in the FISDW phase. 

Interestingly, the presence of the FISDW phase affects 
the Nernst signal at -1 Lebed angle differently. As shown 
in Fig01 at 8 kbar pressure the FISDW onsets (from Rz^) 
coincide the onsets of a large increase of Nernst voltage at 
-1 Lebed angle at various magnetic fields. At still lower 
temperature the voltage at the -1 Lebed angle reaches 
a peak around 300 ~ 400 mK then decreases quickly. 
The peak value is as large as ~ 220 /J.V/K at 7.5 Tesla. 
This is further confirmed by the angle dependence of the 
Nernst effect at base temperature (150 mK), where there 
is large Nernst resonances at ±1 Lebed angles while there 
is none at c'. This behavior is consistent with our previ- 
ous measurements ii£ The effect of FISDW on -1 Lebed 
angles is further confirmed by measurements at higher 
pressures, where FISDW transition temperature Tc(II) 
vary accordingly. For example, at 13 kbar the FISDW 
transition temperature at 7.5 Tesla is suppressed down 
to ^ 350 mK, the enhancement of the Nernst effect at -1 
Lebed angle follows the FISDW transition accordingly as 
shown in the insert of Fig^ 
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The suppression of the Nernst signal at c' by the 
FISDW is not understood at this moment. It is prob- 
ably due to the competition between the FISDW phase 
and "metallic' phase. This difference seems to suggest the 
magic angle c' is different from -1 Lebed magic angle in a 
subtle way. In this paper, we limit our discussion within 
the "normal" state where the MAE is pronounced. We 
note that a full understanding of the MAE should cover 
the FISDW phase, where the MAE is more complicated 
than in the metallic phase. 



D. Field dependence 



tensor S is the product of the resistivity tensor p and the 
thermoelectric tensor a. 
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FIG. 5: Upper panel: Magnetic field dependence of N^x and 
Rzz measured simultaneously at 8 kbar, 1.6 K for c'. Clearly 
Nzx is non-linear with magnetic field just as observed in pre- 
vious measurement^'". 

Lower panel: Ratio of Nza: and Rzz derived from the upper 

Nzx 

panel. As discussed in text, — — oc a^x- The dash line is a 

Rzz 

guide to eyes. It is clear that at low field a^x is linear with 
field. 



As seen in the temperature dependence of N^^; at differ- 
ent magnetic fields, the Nernst effect in (TMTSF)2PF6 is 
very non-linear with magnetic field. In the upper panel of 
Fig|Sl we show simultaneous measurements of magneto- 
resistance R^z and the Nernst effect vs. magnetic 
field at 1.6 K, 8 kbar and 3° off c'. It is clear that the 
Nernst signal has a super-linear field dependence. An 
obvious nonlinear effect is the large magnetoresistance of 
(TMTSF)2PF6. In transport theory, the thermopower 



S = p ■ a 



(1) 




FIG. 6: Magnetic field dependence of N^x and Rzz measured 
simultaneously at 8 kbar and 375 mK for field at 3° away 
from -1 Lebed angle. The upturn of the Rzz around 5.5 Tesla 
defines the threshold field of the FISDW transition, f^zx is 
highly non-linear and enhanced greatly in FISDW phase. 

Therefore, Szx = PzxO^xx ~t- PzyOiyx ~t~ PzzOCzx Pzz^zx- 

Here we ignore the first two terms since the Hall ef- 
fects are negligibly small in the metallic phase for 
(TMTSF)2PF6. To obtain a, we took the ratio of the 
Nernst signal N^^, and the resistance R^^ at the same 

Nzx 

field to obtain field dependence of azx cc — — according 

^zz 

to Eq^ The prefactor depends on the sample geometry. 
The result is shown in the lower panel of FigEl Clearly 
azx is approximately linear with magnetic field below 
5 Tesla. Therefore, the non-linearity of the Nernst sig- 
nal Nz^ in (TMTSF)2PF6 mainly comes from the large 
magneto-resistance. The linear field dependence for field 
along c' and -1 Lebed angle suggests that azx is probably 
a more fundamental quantity in the thermoelectric effect 
of (TMTSF)2PF6. 

FigEl shows the field dependence of Nzx and Rzz for 
field along 3° off the -1 Lebed angle at 375 mK. We can 
see that when the (TMTSF)2PF6 goes into the FISDW, 
the resistance rises up sharply around 5.5 Tesla due to the 
presence of the FISDW gap. The Nernst signal also shows 
a sharp upturn around 5.5 T and rises up dramatically. 
This agrees with the observation of the enhancement of 
the Nernst signal at -1 Lebed angle in the temperature 
dependence (Fig. 



E. The effect of pressure 

The ground state properties of the Bechgaard salts are 
strongly affected by hydrostatic pressure. The tempera- 
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FIG. 7: Zero field temperature dependence of R^z at 8, 10 and 
13 kbar. The T^o's are 1.02 K, 0.87 K and 0.77 K respectively. 



ture, pressure and magnetic field (T-P-H) 3D phase dia- 
gram (Fig.l in RefEOi) summarizes the effects of pressure 
on various phase transitions. For example, the thresh- 
old field of the FISDW phase progressively increases as 
the pressure gets higher. The superconducting transition 
temperature Tc is also slowly suppressed by increasing 
pressure. Figd shows the zero field temperature depen- 
dence of resistance R^z at 8 kbar, 10 kbar and 13 kbar 
respectively. The superconducting transition tempera- 
ture decreases slowly with increasing pressure. Using the 
onset definition(90% of the normal state value), we found 
that the superconducting transition temperatures Tco are 
1.02 K, 0.87 K and 0.77 K respectively In Fig. |Hlwe 

150.0^l 



JOO.OiJ 



> 



50.0|J 



13 kbar/Ztfl'' 



' V \ 8 kbar ■ 
: iTlOkbar 



7.5T 
c'+3° 



0.0 
0.0 



0.5 



1.0 
T(K) 



1.5 



2.0 



FIG. 8: Temperature dependence of the Nernst signal at 
8 kbar , 10 kbar and 13 kbar at 7.5 Tesla. Qualitatively no 
change is observed for various pressures. 

show the comparison of the temperature dependence of 
the normalized Nernst signal at 3° off c' between 8 kbar, 
10 kbar and 13 kbar. Qualitatively, the temperature de- 
pendence of the Nernst signal is pressure insensitive. 
The angular dependence doesn't change significantly 
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FIG. 9: Angle dependence of N^a^C'*") and R2z(thin line) at 
(a) 8 kbar, 2. IK and 7.5 Tesla; (b) 10 kbar, 66QmK and 7.5 
Tesla; (c) 13 kbar, IK and 6 Tesla. Qualitatively, no change 
is observed. 



either as we vary the pressure. Fig. |51 shows the angular 
dependence of N^^; and Rzz at 8, 10 and 13 kbar respec- 
tively. Note that the temperature and magnetic fields of 
these data are not identical. However, qualitatively they 
are all the same. Note that all data shown here are in 
the metallic phase. 



III. BOLTZMANN CALCULATION OF NERNST 
EFFECT 

To gain some elementary intuition about transport 
processes it is often instructive to look first at a general- 
ized Drude approximation, by which we mean a classical 
gas of charged particle in lowest order response to an 
applied set of driving fields. 



A. Drude Transport 

In Fig^l we show a cartoon of the particle motion 
of such a charged gas. In a Drude model forces accel- 
erate particles which then lose momentum in collisions 
at a rate l/r. The basic equation of motion is there- 
fore mv/r = F, the charge per particle is q, the particle 
density is n, the current density is simply charge density 

times velocity j = ngv and j = cr • E. In FigEII (a) 7 F 

9 

g^T" TIC[ T 

is gE and v = — E, j = • E, the conductivity a 

m m 
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is given by cr = . In the presence of the magnetic 
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FIG. 10: Generalized Drude model schematic. Magnetic field perpendicular to the plane, (a), (c) E field drive, (b), (d) VT 
drive, (a),(b), one carrier, (c),(d) two oppositely charged carriers. 



field there is a Lorentz force which deflects particles in 
the y direction building up charges on the upper and 
lower boundaries. The charges continue accumulating 
until the electric field they generate exactly cancels the 
Lorentz force, Ey = Vx ■ Bz- In steady state the Hall 
field Ey — ■ Bz completely compensates the effect of 
the magnetic field, the carriers only drift in the x direc- 
tion and there is no magnetoresistance. In (b) the drive 
is a temperature gradient. In this simple model the gas 
is ideal and we use the ideal gas law, P = nksT. A 
temperature gradient translates to a pressure gradient 

VP 

VP = nksS/T, or a force per particle of F = . The 

n ^ 

charged particles will flow in the x direction charging the 
boundaries and creating an opposing field E^- The cur- 

rent and charging stop when qEx = or E^ = — V^^T. 

q 

This thermoelectric voltage is the Seebeck effect with co- 

efficient Si = — . (This was a big failure of the Drude 
q . 

model. It overestimates the thermopower by several or- 
ders of magnitude. The reason is quite evident today. 
We have a degenerate electron gas (DEG) rather than 
an ideal classical gas. The effective number of degrees of 
freedom, or particles that can transport heat, is reduced 

by ~ — so the Sideg —)• In this picture the 



electric force cancels the pressure and the particles have 
no velocity, v x B = and the Nernst voltage is zero. 
For the simplest conductors we therefore expect the ther- 
mopower to be sizeable and the Nernst effect negligible. 
But it is worth noting that effectively the same argument 
would suggest that the magnetoresistance is negligible. 

In Fig^] (c) and (d) we consider the case of two op- 
positely charged carriers, which are otherwise identical. 

With an electric field along x the two carriers move at 
qr 

Vx — — Ex in opposite directions both contributing to 
m ^ 

the electrical current and conductivity which remains 



nq^T 
m 



Now however, in the presence of B^ both are 

deflected in the same direction, there is no charge accu- 
mulated on the boundaries, there is no Hall voltage, and 



velocities persist in both directions, Vy = v. 
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There is now magne- 
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toresistance, — (uJct) . In (d) the drive is a tempera- 
ture gradient again producing a pressure gradient. Both 
types of particles move down the pressure gradient with 
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velocity, Vx 



m 



there is no charge accumulation, 



and no field generated along x so the Seebeck coefficient 
is zero. In the presence of a magnetic field, the particles 
with the same velocity but opposite sign are separated, 
charges accumulate on the upper and lower boundaries 
until the electric field compensates the Lorentz force, 
xTtBz 

Ey = VxB^ = 



The result is a Nernst volt- 



■ ^. c ks \q\B^T 

age with coefficient byx — = -r—. = biOJcT. 

m \q\ m 

With the degenerate electron gas correction we should 

+1, t- a ksT 

then expect bxyOEG ^ TT~F^y^cT). 

W\ Ef 







Hall 


MR 


Seebeck 


Nernst 


one carrier 


nq^T 


1 

nq 





fcs ksT 





m 


q Ef 


ambipolar 


nq'^T 





Kr)2 





ks ksT 


m 


-r-. ^cT 

\q\ Ef 



The extension of these results to field orientation and 
different ratios of the densities of the oppositely charged 

particles is straightforward. With 5*1 = 



/(«) 



1 + {uj.tY 



\q\ Ef 
1 



a 



we find Si{a) = Sif{a), Sxyia) = Si{ujcT)g{a), Rh = 

-^f{a), — = {LOcTfg{a). The HaU and Nernst volt- 
n\q\ p 

ages vary as B x E and B x VT respectively, /(a) and 
g{a) are plotted in Fig llll Putting in realistic param- 
eters, we will obtain the Nernst effect in Drude model 
is order of 10 nV/K, linear with magnetic field. Drude 
picture only predicts that the Nernst effect has a simple 
sin 9 dependence of magnetic field orientation and linear 
with magnetic field and temperature. 



electron 

1.0 
0.5 ^ 



ambipolar 



hole 



-0.5 



-1.0 







f{a) / 


CO X - 1 



-1.0 -0.5 0.0 0.5 

a={n*-n")/(n*+n") 



1.0 



71 — Tl 

FIG. 11: Drude picture in general. Here a = — ; — ■ is the 



net charge density of mobile charge carriers. 
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B. Boltzmann Transport and QID Fermi Surface 

Simple Drude calculations are not capable of handling 
the highly anisotropic nature of the Bechgaard salts nor 
the angular orientation of the field relative to the lat- 
tice vectors. The simplest treatment which includes 
the bandstructure comes from a steady state Boltzmann 
equation. It has previously been shown that Boltz- 
mann transport, appropriately modified to follow elec- 
tron trajectories over many Brillouin zones (or equiva- 
lently many Umklapp scatterings) can give magic angle 
effect in R^^Moreover, this model is in qualitative agree- 
ment with the measurements in the (TMTSF)2C104 salt. 
Although, Boltzmann transport and the Osada model 
have not been successful for (TMTSF)2PF6, they are 
still the only reasonable single particle treatment avail- 
able. We therefore performed both numerical and ana- 
lytic calculations using a Boltzmann Transport formula- 
tion based on the Tight Binding Approximation band 
structure within the single Relaxation Time Approxi- 
mation. For simplicity,the triclinic crystal structure of 
(TMTSF)2X is taken as orthorhombic. 



-2tn COS k^-a 



2ti, cos kyb 



2tr cos k^ 



(2) 



For (TMTSF)2X, > 4 > U. 
the kx dispersion for simplicity, 
so-called linearized dispersion: 



Often people linearize 
One would obtain the 



e — £/ ~ zLhvf{kx T kf) — 2tbCos{kyb) — 2tcCos{kzc) (3) 

In this approximation, the Fermi velocity vj (or the den- 
sity of states on Fermi surface Nf{ef) = 1/hvf) is a con- 
stant for a given energy. Many of the angular magneto- 
resistance oscillations (AMRO), e.g. the Danner-Kang- 
Chaikin oscillation (ac-rotation)^, the third angle effect 
(ab-rotation)^ and the combination of themS can be 
understood within this approximation. There is excellent 
agreements between experiment and theory, especially 
for (TMTSF)2C104. However, for a linearized dispersion 
relation the Hall effect is 7.eio{axy = <Jxz = 0)- It is not 
surprising the Nernst effect is also zero {Sxy = Sxz=0) 
in this approximation. In order to calculate the Nernst 
effect, we have to use the full dispersion or a nonlinear 
approximation in either numerical computation or in an 
analytic calculation. 

In general transport theory, we consider both an elec- 
tric current J and a thermal current Jq in response to 
an electric field E, a magnetic field B and a temperature 
gradient (-VT): 



J = cr-E + a- (-VT) 
Jg = Ta' • E + K- (-VT) 



(4) 



Here a is the electric conductivity tensor, a and a' the 
thermoelectric tensor and k is the thermal conductiv- 
ity tensor. Here Q;ij(H) = a^-j(— H) according to On- 
sager relation. ^'^'^^ For free electron gas, these coeffi- 
cients can be obtained by applying the relaxation time 
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approximation 



,36 



Numerical Calculation 



47r3 



dSk 

hv 



I 7r\ 



v(v) 
2 da{e) 



de 



(5) 



Here we assume an energy and momentum independent 
relaxation time r. The velocity average is defined as: 
dt 



(v(k(t))) 



v(k(t)), where k(i) is the semi- 



classical motion of electrons on Fermi Surface in the pres- 
ence of a magnetic field. Eqs. [51 are the starting point 
of our calculation. From the general transport equations 
(Eqs. 21, by setting J = as the boundary condition in 
thermoelectric measurements, one would find EqQ The 
Nernst signal N^^; is an off-diagonal element of thermo- 
electric power tensor S. 

Now let's focus on a quasi-lD (QID) system, 
(TMTSF)2X. The Fermi surface of QID consists of two 
slightly warped sheets. In order to evaluate the Fermi 
surface averaging velocity (v) we have to calculate the 
motion of electrons on Fermi surface in the presence of 
the magnetic field. Here the a-axis {kx) is the best con- 
ductivity direction, so the Fermi surface is approximately 
normal to the k X Exis. TlicreforG, hx is 9- function of ky^ 
kz and e from the dispersion relation (Eq. [21), i.e. 



kx — kx {ky, kzf s) 



(6) 



Therefore, kx is not an independent variable for the semi- 
classical motion of electrons, given that the motion of the 
election is confined to the Fermi surface in the presence 
of the magnetic field B. The Equations of motion can be 
reduced to: 



dt 

dkz_ 
dt 



VyBx 



(7) 



By solving these two equations of motion EqQ we can 
evaluate v(t) = v(k(i)). This is what we need for the 
conductivity tensor a. ^From Eqs. [Sj the thermoelectric 
coefficient tensor a is proportional to the energy deriva- 
tive of the conductivity tensor (7(e) at the Fermi energy. 



dky dk^ 



h\vx 



klerT d 
127r de 



(8) 



dky dk z 



h\vx 



Vi (Vj) 



If we could find an analytic form of cr{s), evaluating a 
would be straightforward. In our analytic calculation, 
we derive an approximate analytic form for cr{e), and we 
obtain a by taking differentiation. On the other hand, it 
is straightforward to evaluate both a and a numerically. 
Here a can be calculated by taking the energy derivative 
of each term in the integral. (EqlSJ 



From Eq. |Slwe can obtain: 
klerT , , 

dkydk^ 



de V \vx\ J \vx\ de \vx\ de 



(a) 



(b) 



Note here a cx T if we assume r is T independence, 
which is a good approximation at low temperature. It is 
straightforward to evaluate the first two terms. Since: 



— - h M^^ — 
de de 

Here we define the inverse mass tensor A4~^ as: 

ldui_ _ 1 d^ejk) 
h dki h? dkidki 



(9) 



(10) 



Here ky and kz are independent variables, so the energy 
dependence of the first two terms only comes from kx, 
which is a function of e (see Eql^J. Then, we will have: 



For 





d 


dkx 


d 


1 d 




de 


'dT 


dkx 


hvx dkx 


d 

Therefore: — 

de 






7,3 


dv, m~l 
and — = S^^x^^ 
de Vx 



the (c) term, we exchange the differentiation and aver 
d (v) / dv 



aging(integral), i.e.- 



de 



— ). Putting the (a), (b) 
de ' 



and (c) terms together, we obtain an expression for the 
thermoelectric coefficient tensor a: 



a — a 

where a'^^-' is defined as 
fc2 erT 



(+) 



(12) 



a 



(±) 



127rfi. 



dkydkz{i) 



dv 



dv^ 
'dF 



kx<0 



Now we need to find — — (t) , which mvolves the motion of 
de 

dv 

an electron in the magnetic field, ^^(i) depends on the 

de 

value of the Fermi Energy not only through kx , but also 
through ky and kz, because ky and kz are also functions 
of t when the electron is moving on the Fermi surface in 
the magnetic field. These functions depend on energy e 
and the "initial" condition (fc°, /c°), i.e. 



ky{t) — ky{t; e, fc°, /c^) 
kz{t) = kz{t; e, fcJJ, fc^) 



(13) 
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FIG. 12: Angular dependence ol conductivity tensor a. Here B=8 Tesla. The origin is the c-axis, while 90°is b-axis and 
30°is the +1 Lebed angle: c+b. The inset of ct^z shows a hump at the +1 Lebed angle. The open circles are the numerical 
calculations, and the thin hues are the analytic results. 



We can't derive the energy dependence of ky and 
directly if we don't know the solution of the Equations of 
motion. However, we can find the differential equations 

which are satisfied by and -r-^, from the equations 

OS oe 
of motion of ky and fc^ (EqEIl. 



d dky 
dt de 

d dky 
dt de 



d dky 
de dt 

d dkz 
de dt 



dvx 
de 



dvz 
de 

de ^ de ^ 



(14) 



From (Eq|SJ), we will find 

dk,{t) _ 1 



de 



hvx{t) 



1 — hv,, 



dky it) 
de 



hv 



dkS) 
de 



Combining Ea ll4l with Eq|^ wc can obtain numerical 
d'v 

solutions of -^{t). Since ky, k^ are independent of e at 



t = 0, the "initial" conditions for -:r^(t) and -^(t) are 

de de 

In summary, to calculate both a and a, we numerically 
solve two sets of equations of motions, Eqs. [7|and Eqs. 
n~n To treat the differential equations in Eqs. {T] and 
[T^ . we use a 4th order Runge-Kutta Method.^''' Then 

' dv 



we numerically integrate (v) and \^'q~ J ^^^^ evaluate 

both a and a by Fermi surface integrals. Here we use 
the band parameters {ta — 0.25eV, tb = O.Q2AeV and 
tc — O.OOSeV) from tight binding approximations, realis- 
tic lattice parameters (a=3.49 A, h—7.7 A and c=13.264 
A) and a scattering time r = 4.26 x lO^^^sec from previ- 
ous studies by Banner et alm^. Here we use B=8 Tesla, 
and T=l K, which are comparable with experiment con- 
ditions. The combination of the parameters produce 
cBt 

LOcT = = 6 for an isotropic free electron gas. We use 

me 

a 20 X 20 grid on the Fermi surface in the calculations. 




Although it is a httlc rough, it catches aU the main fea- 
tures. We also used a 40 x 40 grid for some points and 
did not find a significant difference. 



D. Analytic Calculation 



The basis of our analytic approximation scheme is find- 
ing the proper correction to the linear dispersion approx- 
imation Eq|2| We expanded to next order to include 
the effect of non-linearity of the dispersion, i.e. ~ + 
Sv / . This approach is basically the semi-classical version 
of Lebed's quantum approach'^^. It is straightforward to 

find: Vf — — — smfe/fl 



Svf_ 



h 

cos kfa 
sin^ kfa 



here coskfa = —-r— and 



tb 



2ta 



cos{kyb) + — cos{kzc) 



„ „ . „ cosfc/a tb 
Definmg (3 = — . „ — 
sm kfttta 



and 7 = 



cos kfa t, 
sin^ kfat, 



then 



we have: 



Vx ^ Vf[l ~ f3cos{kyb) — ^cos{kzc)] 
2tbb . 

jy = sm kyO 



(15) 



2tcC 



sin fcrC 



By substituting v^. in the equations of motion Eqs|3 
one can obtain analytic expressions of ky{t) and k^it) 
and evaluate velocity averages (v) for <j{e). Then it is 
straightforward to obtain a from Eq|Sl Details of the 
analytic calculation are presented in the appendix. As 
shown in Fig^l Fig^l FigtTH and Fig^l our analytic 
calculation reproduces the main features of the numeri- 
cal calculation. In some cases, the results from different 
methods overlap. Therefore, we believe our calculations 
describe the main behavior of the Nernst signal in Boltz- 
mann transport within the tight binding approximation. 
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FIG. 14: Thermoelectric coefRcient tensor a. In Boltzmann transport, a ~ ^ 
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E. Results 

FigEl shows the calculated angular dependence of 
the conductivity tensor a. The graphs are arranged 
in the pattern of the tensor elements aij in the matrix 
form. Here 6 is defined respect to the c-axis. There- 
fore 9 =90°corresponds to the magnetic aligning at the 
b-axis. As shown in the insert graph of azz, there is a 
small hump at 1st Lebed angle: c+b, which is about 
30°in our approximation. This is one test that our calcu- 
lations reproduce the angular dependence of a^z calcu- 
lated by Osada et alS^. By increasing the scattering time 
T, or magnetic field, we can clearly resolve a peak at this 
angle. We also confirmed other AMROs, e.g. ac-rotation 
and reproduce the Banner oscillations^!. 

However, as far as we know, there is no calculation 
of tensor a in the literature for comparison. This was 
our motivation for performing the analytic calculation to 
confirm our numerical results. Once we obtained a and 
a, we got the thermoelectric power tensor S fFig ll5|l by 
taking the product oi p = a^^ (Fig. I13|) and a (Fig. I14|l . 



It is clear that the quality of numerical calculation of 
a is much better than that of a. For a most curves are 
very smooth and only minor oscillations are observed. 
For a most curves are smooth, except Uyz and azy oiyz 
has some spiky features for 9 close to c; while azy has 
some spiky features for 9 close to b. This is because 
the energy derivative of the velocities is very sensitive to 
the location on the Fermi surface (fcj,, kz). Finite size 
grid integration could also generate artificial spikes if the 
integrand oscillates a lot. 

By varying the grid size and the integration cut off 
limit, these artificial features can be suppressed, but at 
the expense of much more computation time. Since we 
are only interested in the general behaviors and magni- 
tudes for a given set of parameters, we will use these 
non-perfect calculation results to compare with experi- 
ments, while keeping in mind that sharp features might 
be artificial. Also, our analytic results will help us to find 
out the physical features. 

The Nernst signal Szx corresponds to the experimental 
results discussed in the previous section. It is clear that 
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Thermopower Tensor S 




FIG. 15: Thermoelectric Power tensor S, which is the product of resistivity tensor ofFig lT^ and thermoelectric coefficient 
tensor a (Fig^J. 

its angular dependence is sin^-like, which agrees with 
the simple Drude model. The maximum magnitude is 
about 1 fiY/K, which is 2 orders smaller than what we 
found in (TMTSFjaPFg as shown in Fig. [TBI This result 
is very different in shape from our observation, missing 
the resonances at magic angles and it gives the wrong 
temperature dependence (5 oc T in Boltzmann Trans- 
port). Therefore, we conclude that the single particle 
picture is not able to describe the Nernst effect observed 
in (TMTSF)2PF6. 

It is interesting to note that the Nernst signal Sxz has 
a similar angular dependence as the experimental Nernst 
effect Szx , showing a peak near c and a sign change. How- 
ever, the geometry is completely opposite and the value 
is about 8 orders of magnitude too small. Experimen- 
tally we couldn't detect a sizable Sxz, though we did not 
optimize the experimental setup for that measurement. 
Comparing the off-diagonal elements of 5 {i.e. the Nernst 
effect), we find that the elements in the upper triangle, 
Sxy, Sxz and Syz are much smaller that those in the lower 
triangle, Syx, Szx and Szy. More surprisingly, the conju- 
gate elements don't have the same angular dependence. 



Does this violate the Onsager relation^^iSi? The answer 
is no. The Onsager relation only states constraint on a 
and a and a' (defined in Eqs. 

aij{H) = aji{-H) , . 

a,,(/f)=a;,(-i7) ^^^) 

By symmetry, we know that aij (aij) is an odd func- 
tion of magnetic field H for i ^ j. 

= -cr^j{-H) , , 

Together with the Onsager relation Ea ll6l we can find: 



cfij{H) = -aj,{H) 

This is exactly what we see in the calculations (except 
ayz (ctyz) and fJzyioizy), which are not real Hall effects 
since the magnetic field is in the plane). However, the 
thermopower tensor S is the product of p = a^^ and a. 
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FIG. 16: Comparison of Nz^^between calculations (Fig. I15|l 
and experimental data (Fig. |5J. 



IV. DISCUSSION 

The giant value of the Nernst effect and the Nernst 
resonances at Magic angles are not understood and it 
appears difficult to explain them in conventional Fermi 
Liquid models as illustrated by comparing our experi- 
ments and Boltzmann transport calculations. An exotic 
feature is that the Nernst signal changes its sign sharply 
at magic angles, with 3 "resonances" within 70° in our 
measurements. As we know, the sign of the transverse 
electric field is determined by the cross product: VT x B. 
Of course, the physics really involves E = v x B. Since 
the temperature gradient is fixed, the only quantity that 
could possibly change its sign is the component of B. As 
the magnetic field passes through a magic angle, the only 
component of magnetic field that could change sign rela- 
tive to a magic angle is the one that is perpendicular to 
the direction of the magic angle. Therefore, we have to 
conclude that the Nernst signal in (TMTSF)2PF6 comes 
from V X B^. This means that whatever is moving is 
confined in the Magic angle plane. The main idea that 
underlies our interpretation of the Nernst resonances is 
that the transport is only coherent at the planes which 
are "parallel" (or close to parallel) to the magnetic field. 
In other words, the coherent electronic motion is con- 
trolled by the orientation of the magnetic field relative 
to the planes defined by the conducting chains and the 
interchain directions. (Fig ll7|l 



In general one should not expect Sij = —Sji. This is 
only true when we consider an isotropic system, where 
Pa = Po and pij = pn for i =/= j. In an anisotropic 
system like (TMTSF)2X, pr^^ ^ Pyy ^ Pzz- If we ignore 
the Hall effect, we will have Sxy ~ PxxOtxy ^ Pyy^yx ~ 
Syx- Nevertheless, the Nernst effect Sij {i ^ j)is an odd 
function of H and Seebeck effect Su is an even function 
of H as long as there is a inversion symmetry. 



In summary, we numerically and analytically calcu- 
lated the thermopower tensor iS by evaluating both the 
conductivity tensor a and thermoelectric coefficient ten- 
sor a. The numerical results agree well with the an- 
alytic approximation. This gives us confidence on the 
reliability of our calculations. It is clear that Boltzmann 
transport within a single particle picture is not consistent 
with our observation in (TMTSF)2PF6. Therefore, cor- 
relation effects due to the strong e-e interaction should 
be considered in understanding the giant Nernst effect 
found in (TMTSF)2PF6. However, we note that the an- 
gular dependence of S^x fits the data in (TMTSF)2C104 
very well, though there is a factor of 10 or so difference 
in magnitude's Our results are not limited to the Bech- 
gaard salts (TMTSF)2X. For any QID system with open 
Fermi surface, all the transport coefficients can be calcu- 
lated using our results based on Boltzmann transport in 
a tight binding model. Our original results should prove 
useful for further investigations. 




FIG. 17: The schemaic picture indicates Coherent Transport 
only in the magic angle planes which are parallel to magnetic 
field. The Nernst signal would chagne sign as the field is 
rotated above or below these planes. | -1 | : B |[ -IL, a(c'-b')- 
plane; | c' | : B || c', ac'-plane; | -1-1 | : B || -I-IL, a(c'-|-b')-plane. 
The Nernst data is taken from previous measurements—. 

The nature of the coherence is not clear at this 
moment. One possibility is quasi-particle coherence. 
For example we may have a Field Induced Inter- 
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chain/ Interplane Decoupling picture. 



CD 





FIG. 18: Field Induced Interchain /Interplane Decoupling 

[T] : Bx > B*, = 0, chains are decoupled, Nernst signal 
is zero ejv = 0; 

\k 

0: 



Bx < B*, tc'^ / 0, chains are coupled, ejv > 0; 
Bx = 0, ^ 0, chains are coupled, ejv = 0; 
Bx < B*, tc'^ 7^ 0, chains are coupled, ejv < 0; 



Bx > B*, i 
is zero ejv = 0. 



eff 



0, chains are decoupled, Nernst signal 



FigEl shows the basic idea of this picture. 
(TMTSF)2PF6 is a QID system, which consists of con- 
duction chains along the a-axis. The chains are weakly 
coupled with each other. When the magnetic field B 
is far from a magic angle, e.g. c', the inter-chain coher- 
ent coupling in this direction is effectively suppressed by 
the large normal field component, i.e., = when 
Bj^ > B* (here B* is a cross-over magnetic field scale). 
In other words, there is no coherent transport in the ac- 
plane for sufficient perpendicular field. (FigE|[T], [5] 
) When B is parallel to the c' direction, the coherent 
coupling in the c-direction is restored and a coherent ac- 
plane is formed. (Fig^l 3 ) If we tilt the magnetic field 



slightly upward, there is a small component B_l of B nor- 
mal to this plane. (FigllSI 2 ) Here is small enough 



that B^ < B*. So the plane is still coherent. In the pres- 
ence of temperature gradient — VT, there is a transverse 
electric field E, the Nernst effect. When we tilt the field 
slightly to the other side of the magic angle, everything is 
the same except the sign of B^ reverses. Therefore, the 
sign of the Nernst signal is reversed. (Fig^l[T] ) This 
scenario is then repeated at the other magic angles as in 
FigE] 

In fact, the idea of field induced decoupling is not new. 
Strong et alm^ considered the isolated conducting planes, 
ab-planes, of (TMTSF)2PF6 as 2D Non-Fermi-Liquid 
due to the strong e-e interactions. The possible non 
Fermi Liquid ground state of "isolated" (TMTSF)2PF6 
chains is supported by transport, optical^ and thermal 
transport measurements. The effect of field induced co- 



herent coupling/decoupling is supported by temperature 
and angular dependence magneto-resistance studies'*—. 
However, this theory is not universally accepted due to 
the unknown nature of the non Fermi Liquid state and 
the lack of a detailed model. 

Another possibility is the 2 dimensional superconduct- 
ing phase coherence proposed by Ong et al.^'^ . Ong 
points out that normal quasi-particles usually give a ther- 
mopower much larger than the Nernst signalfSfl whereas 
we have a large Nernst signal with undetectably small 
thermopower. On the contrary, vortex flow naturally pro- 
duces an electric field that is predominantly transverse. 
This is generally true for most conventional systems, as 
well as HTc cuprateSfSfl Implicit in this model is the abil- 
ity of the magnetic field to destroy phase coherence in 
the planes to which it is normal. It is natural that the 
vortices penetrating perpendicular planes destroy super- 
conductivity. This is similar to the decoupling model 
discussed above and has the consequent" resonances" at 
magic angles. 



V. CONCLUSION 

We present a detailed study of the Nernst effect N^^, 
in (TMTSF)2PF6 as a function of temperature, pressure 
and magnetic field magnitude and direction. The data 
agree well with our previous measurements 4S We have 
calculated Boltzmann transport coefficients by both nu- 
merical and analytic methods with realistic band pa- 
rameters within the single relaxation time approxima- 
tion. The two calculations agree with either very well, 
but fail to describe the experimental data. The large 
magnitude, resonant-like angular field dependence and 
the non-linear field and temperature dependence cannot 
be understood within the semi-classical approximation of 
Boltzmann transport. The sign change of the Nernst ef- 
fect at magic angles strongly suggests that the transport 
is effectively 2-dimensional in lattice planes parallel or 
close to the orientation of magnetic field. The nature of 
the coherence is not clear at this moment. We suggest 
that the answer may lie in field induced decoupling of the 
strongly correlated chains. 
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VII. APPENDIX 



Define (pb — kyb and (pc = kzC In the presence of a 
magnetic field in the bc-plane B = (0, Bsin^, i?cos0). 



define 



ujb cos 9 and 



,s — 



LOr sin 9, here huih 



\e\vfBh and hujc = \e\vfBc. e < is the electron charge. 
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Then the equation of motion become: 



dt 



= uj^{l - (3cos(t>h ~ 7C0S (t)c) 
— — t^c(l ~ /3cos0fc — 7 cos 0c) 



(19) 



The exact solutions are hard to obtain. However, since 
/3, 7 ^ 1, the solutions can be approximated to 1st order 
by an iterative method: 

cj,,(t) ^4>l+ ult ~ /3sin(0° + u^lt) - 7^ sin(0° - J,t) 

Mt) = 0c - ^ct + 13^ sin(0SJ + utt) - 7sin(0° - u^.t) 

(20) 

Note here — |- diverges as 6* ^ 0, and — §■ diverges as 

6* -> 90°. So this solution is only good at < 6^ < 90°. 
One has to be aware that as B approaches b or c, this 



solution may not give the correct result. With the help 
of Jacobi's expansions, it is straightforward to evaluate 
v(t) fEa ll5|) and (v). Then we can obtain analytic ex- 
pression for a. After tedious but straightforward calcu- 
lation, we obtain the conductivity tensor tr. Here we 
assume 7^/3^1, and only keep the leading terms in 
(3, 7, 0^ and (3^. For simplicity, we also use the anti- 
symmetric property of a. It is clear from Fig ll2l that 
the analytic calculation reproduces the numerical results 
very well. Due to the limitation of our expansion, not 
every minor detail was reproduced. For example, in this 
1st order approximation, Gxx is independence of angle 
in analytic form, while numerically it shows a very weak 
angle dependence. The good agreement between differ- 
ent calculations give us confidence about the reliability of 
the calculations. Once we know the analytic form of the 
conductivity tensor cr, we can take its energy derivative 
to obtain a. 



J 



^^(l-/3Ji(/3)-7Ji(7)) 

TTft be 



7Th?C t 1 + (w^t)2 

i^b - ^c)r (w. 



1 + (2c.«r)2 



1 {luI - w«)2t2 1 + {iol + )2r2 



'xy 



'yx 



Ae^Ttl b 



TTh^Vf c\l + (uJ^t)^ 



2 1 + (2c^«t)2 



1 + {uj^T - ^Irf 1 + (i^lr + c^«t)2 



Y^^(^^ + (2c.c^r)2 



7 i^^ 



2_2 



=2, 



TTh-^Vf \ U!^ 

— ^3 ^'^1 "f 6 



1 + {lu^^t + Lu^^ry 1 + (c^»r - c^eT)2 



1 + (u;«r - u;eT)2 1 + (c^^r + c^»t)2 
iJiil) 1 



'ie'^Ttl c 

irh^Vf b\l + (wSr)2 2 1 + (2w«t)2 

6 r 



1 + {lu"^t - a;«r)2 1 + (cc;«r + lo^^t)^ 



(21) 
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^ {(1 - 13 Jm - iMl)) + 3[/3Ji(/3) + 13^ Jm + iMl) + l^J[{l)]} 



m'^c tl [1 + (a;^T)2]2 



a. 



27rA;?,TeT r ,a;«r[l + 2(u;fT)2] 73 



[l + (c.«r)2] 



4 1 + {2ujItY 



' , , 5 + 7(2w^r)2 ^ 



1 + (2a;^) 



7^ 
2 ' 



/3Ji 7-1 +7^1 /?-! 



+ ■ 



37/3 

4 



|l + K'T + .^»r) 



212 



+ 



l + (a;^-u;^)2 ' 1 + (c.^^ + c.^)^ 



-a. 



27rfc|rer&tg l + 3(u;gT)^ 
3/i2oc ^2^1^ (a;^r)2]2 



TTk%TeTl3j'^ 
3^2^^ 

'2a;,^ 



^1 



2 + 3KV - uj'.t)^ 2 + 3(a;f r + u^^rf 



+ 



[l + (a;,^-a;,^T)2]2 [1 + (a;,^ + c.M2 
1 



+ 



1 + (wf T - a;,^r)2 1 + (a;f r + w^r)2 



[l + K^-a;^)2]2 [l + (a;,^ + a;fr)2]2 



+ 



Oizz = 



¥ab 



1 + (w^T - w^t)2 ' l + (a;fr + a;,^r)2 



f l + 3(c^gr)^ 7J1 
\[l + (a;^)2]2 2 



7Ji(7) 7+ll.(2a;^)2 37V((7) 



2 [1 + (2a;^r)2]2 1 + (2a;^T)2 
2 + 3K«r-a;^T)2 2 + 3(a;f r + r)^ 



[1 + {wtr - a;,^r)2]2 [1 + {oj^t + a;^)2] 



(22) 



+3/3^ J( 



1 



1 + (wj^ - r)2 l + (a;^V + a;^r)2 



Once we have a and a, we can calculate iS = • a. 
The analytic form of <S would be too long to write down 

here since every term in S is the sum of three products 
of two matrix elements. We just numerically calculate 



the matrix product. Here we only present the angular 
dependence of all quantities in order to compare with 
experiment. In principle one could obtain temperature 
and field dependence with these formulae. 
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